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Summary. Here we report the discovery of very high energy (> lOOGeV) 7-ray 
emission from the radio emitting X-ray binary LS I +61 303 with the Major Atmo- 
spheric Gamma Imaging Cherenkov (MAGIC) telescope. This high energy emission 
has been found to be variable, detected 4 days after the periastron passage and 
. lasting for several days. The data have been taken along different orbital cycles, 

' and the fact that the detections occur at similar orbital phases, suggests that the 

£SJ , emission is periodic. Two different scenarios have been involved to explain this high 

0^ ' energy emission: the microquasar scenario where the 7-rays are produced in a radio- 

U ' emitting jet; or the pulsar binary scenario, where they are produced in the shock 

which is generated by the interaction of a pulsar wind and the wind of the massive 

■ companion. 

1 Introduction: The MAGIC Telescope 

O , 

The Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescope U is a very 
C/3 . high energy (VHE) 7-ray telescope, operating in a energy band from 100 GeV to 10 

TeV, exploiting the Imaging Air Cherenkov (IAC) technique. Located on the Canary 
Island of La Palma, at 28°45'30"N, 17° 52' 48"W and 2250 m above sea level. The 
telescope has a 17-m diameter tessellated parabolic mirror, and is equipped with a 
3.5°-3.8° field of view camera. See |5j for a complete description of the instrument. 



2 The 7-ray binary LS I +61 303 



LS I +61 303 belongs, together with LS 5039 |3] and PSR B1259-63 g], to a new 
class of objects, the so-called 7-ray binary systems, whose electromagnetic emission 
extends up to the TeV domain. In section 12.31 we will describe the possible VHE 
emission scenarios for this object: it was proposed as microquasar candidates, but 
also thought to be in a pulsar binary scenario. Thus, the question of whether the 
three known 7-ray binaries produce TeV emission by the same mechanism, and by 
which one, is currently object of an intense debate [5J. 
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2.1 Description of the object 

This 7-ray binary system is composed of a BO main sequence star with a circumstellar 
disc, i.e. a Be star, located at a distance of ~2 kpc. A compact object of unknown 
nature (neutron star or black hole) is orbiting around it, in a highly eccentric (e = 
0.72 ±0.15) orbit. 

The orbital period -with associated radio JBJ and X-ray [7J outbursts- is 26.496 
days and periastron passage is at phase 0.23 ± 0.02 

Radio outbursts are observed every orbital cycle at phases varying between 0.45 
and 0.95 with a 4.6 years modulation 

High-resolution radio imaging techniques have shown extended, radio-emitting 
structures with angular extension of ~0.01 to ~0.1 arc-sec, interpreted within the 
framework of the microquasar scenario, where the radio emission is originated in a 
two-sided, possibly precessing, relativistic jet (/3/c = 0.6) |10|. However, the two- 
sided jets were not completely resolved and no solid evidence of the presence of an 
accretion disk (i.e. a thermal X-ray component) has been observed. There are hints 
of variability of the 7-ray flux 

LS I +61 303 was considered as one of the two microquasar candidates posi- 
tionally coincident with EGRET 7-ray sources |12|. and the only one located in the 
Northern Hemisphere -hence a suitable target for MAGIC. The large uncertainty of 
the position of the EGRET source did not allow an unambiguous association. 

2.2 MAGIC observations 

LS I +61 303 was observed with MAGIC for 54 hours (after standard quality se- 
lection, discarding bad weather data) between October 2005 and March 2006 |13| . 
As mentioned above, MAGIC is able to operate in moderate moon conditions, and 
in particular, 22% of the data used in this analysis were recorded under moonlight. 
The data analysis was carried out using the standard MAGIC reconstruction and 
analysis software |14l 1151 ITC) . 

The reconstructed 7-ray map is shown in Figure The data were first divided 
into two different samples, around periastron passage (0.2-0.3) and at higher (0.4-0.7) 
orbital phases. No significant excess in the number of 7-ray events is detected around 
periastron passage, whereas there is a clear detection (9.4<r statistical significance) 
at later orbital phases. 

The distribution of 7-ray excess is consistent with a point-like source located at 
(J2000): a = 2 h 40 m 34 s , 5 = 61°15'25", with statistical and systematic uncertainties 
of ±0.4' and ±2', respectively. This position is in agreement with the position of 
LS I +61 303. In the natural case in which the VHE emission is produced by the 
same object detected at EGRET energies, this result identifies a 7-ray source that 
resisted classification during the last three decades. 

Our measurements show that the VHE 7-ray emission from LS I +61 303 is 
variable. The 7-ray flux above 400 GeV coming from the direction of LS I +61 303 
(see Figure has a maximum corresponding to about 16% of the Crab nebula flux, 
and is detected around phase 0.6. The combined statistical significance of the 3 
highest flux measurements is 8.7<r, for an integrated observation time of 4.2 hours. 
The probability for the distribution of measured fluxes to be a statistical fluctuation 
of a constant flux (obtained from a \ 2 fit of a constant function to the entire data 
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Fig. 1. Smoothed maps of 7-ray excess events above 400 GeV around LS I +61 303. 
(A) Observations over 15.5 hours corresponding to data around periastron (phase 
0.2-0.3). (B) Observations over 10.7 hours at orbital phase 0.4-0.7. The position of 
the optical source LSI +61 303 (yellow cross) and the 95% confidence level contours 
for two EGRET sources are shown. From Albert et at |13|. 

sample) is 3 x 10~ 5 . The fact that the detections occur at similar orbital phases 
hints at a periodic nature of the VHE 7-ray emission. 

Contemporaneous radio observations of LS I +61 303 were carried out at 15 GHz 
with the Ryle Telescope covering several orbital periods of the source. The peak of 
the radio outbursts was at phase 0.7, i.e. between 1 and 3 days after the increase 
observed at VHE 7-rays flux (see Figure [IJ. 

The VHE spectrum derived from data between ~200 GeV and ~4 TeV at orbital 
phases between 0.4 and 0.7 is shown in Figure The obtained spectrum is fitted 
well (x 2 /ndf = 6.6/5) by a power law function: 

dN 1 /{dA/dt/dE) = (2.7 ± 0.4 ± 0.8) x 10~ 12 x S (-3-6±o.2±Q.a) cm - 2g -i TeV -i 

where iV 7 is the number of 7-rays reaching the Earth per unit area A, time t and 
energy E (expressed in TeV). Errors quoted are statistical and systematic, respec- 
tively. This spectrum is consistent with that measured by EGRET for a spectral 
break between 10 and 100 GeV. 

We estimate that the flux from LS I +61 303 above 200 GeV corresponds to 
an isotropic luminosity of ~ 7 x 10 33 erg s _1 (assuming a distance to the system 
of 2 kpc |18|L The intrinsic luminosity is of the same order of that of the similar 
object LS 5039 and a factor ~ 2 lower than the previous experimental upper limit 
(< 8.8 x 10 -12 cm" 2 s _1 above 500 GeV) obtained by Whipple [S]. LS I +61 303 
displays more luminosity at GeV than at X-ray energies, a behavior shared also by 
LS 5039. 

2.3 Emission scenarios 

Microquasars are a subclass of stellar, X-ray binary systems that display prominent 
radio emission, usually attributed to the existence of jets of relativistic particles. 
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Fig. 2. VHE 7-ray flux of LS I +61 303 as a function of orbital phase for the 
six observed orbital cycles (six upper panels, one point per observation night) and 
averaged for the entire observation time (bottom panel). Vertical error bars include 
lcr statistical error and 10% systematic uncertainty on day-to-day relative fluxes. 
Only data points with more than 2a significance are shown, and 2a upper limits |17| 
are derived for the rest. The modified Julian date (MJD) corresponding to orbital 
phase is indicated for every orbital cycle. From Albert el al. |13|. 



They are named after the similarities with active galactic nuclei (AGNs), since 
microquasars show the same three ingredients that make up radio-loud AGNs: a 
compact object, an accretion disc, and relativistic jets 12(11 . Hence, microquasars are 
galactic, scaled-down versions of an AGN, where instead of a super-massive black 
hole we deal with a compact object of just a few solar masses that accretes material 
from a donor star. The similarities with AGNs explain the large interest risen by mi- 
croquasars. In fact, the short timescale variability displayed by microquasars allows 
to see changes in the ongoing physical processes within typical time scales ranging 
from minutes to months, in contrast with the usual scales of years to observe such 
variability in AGNs. In addition, microquasars could measurably contribute to the 
density of galactic cosmic rays 
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Fig. 3. Left: LS I +61 303 radio flux density at 15 GHz measured with the Ryle Tele- 
scope (green squares, right axis) and results from the last orbital cycle observed by 
MAGIC (black dots, left axis). Right: Differential energy spectrum of LS I +61 303 
for energies between 200 GeV and 4 TeV and averaged for orbital phases between 
0.4 and 0.7, measured by MAGIC. The dashed line corresponds to the Crab nebula 
differential spectrum also measured by MAGIC. The solid line is a fit of a power 
law to the measured points. 



LS I +61 303 was usually thought to be similar to the microquasar LS 5039, be- 
cause evidence of relativistic jets has been found at radio frequencies. In this scenario, 
the high energy emission would be produced in shocks at the relativistic jets |20|. 
The observation of jets has triggered the study of different microquasar-based 7-ray 
emission models, some regarding hadronic mechanisms: relativistic protons in the 
jet interact with non-relativistic stellar wind ions, producing 7-rays via neutral pion 
decay; some regarding leptonic mechanisms: IC scattering of relativistic electrons in 
the jet on stellar and/or synchrotron photons |22j . 

However, no clear signal of the presence of an accretion disk (in particular a 
spectral feature between ~ 10 and ~ 100 keV due to the cut-off of the thermal 
emission) has been observed so far. Because of that, it has been alternatively pro- 
posed that relativistic particles could be injected into the surrounding medium at 
the shock where the wind of the young pulsar and the wind of the stellar companion 
collide |23|. which seems to be the case of PSR B1259-63. 

In the case of LS I +61 303 the resemblance of the time variability and the 
radio/X-ray spectra with those of young pulsars support such hypothesis. However, 
no pulsed emission has been detected from LS I +61 303. 

Observation at VHE with MAGIC simultaneously together with other instru- 
ments at other wavelength domains -in particular radio- will help elucidate the 
mechanism of the TeV emission and hence the nature of 7-ray binaries. 
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